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Abstract 

We calculate numerically the renormalization group (RG) flow of lattice QCD in two- 
coupling space, iPixi, f3ix2)- This is the first expHcit calculation of the RG flow of SU(3) 
gauge theory. From the RG flow, a renormalized trajectory (RT) is revealed. Its behavior 
is consistent with the strong coupling expansion near the high-temperature flxed point. 
Actions with (/3ixi, /3ix2) are studied; the lattice spacing is evaluated by measuring the 
string tension from the heavy quark potential. Recovery of the rotational symmetry is 
studied as a function of the ratio /3ix2//3ixi- 
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1 Introduction 



Since Wilson's first numerical RG analysis of SU(2) gauge theory there have been 
many Monte Carlo RG studies of non-perturbative /3-functions (see Ref . Q and references 
therein). In these analyses indirect information about the P function, such as A/3, has been 
obtained [||. Recent progress of lattice techniques [@, ||, ||] allows us to estimate directly 
the renormalization group (RG) flow in a multi-coupling space 0. Here we will discuss RG 
behavior of lattice QCD in the quenched approximation, i.e., pure SU(3) gauge theory. 

New blocked actions 5" as a function of blocked link variables U^^s are constructed 
from the original S{U) as 

e-^'(^") = I e~^(^)<5(C/^ - PmdU, (1) 

where P defines the blocking transformation. S and S' are different points in the infinite 
coupling space: after the blocking transformation S moves into S". Only fixed points are 
left invariant. In this process the lattice spacing, i.e., the cut-off of the theory, is changed, 
and this RG flow characterizes the theory. 

Corresponding to each blocking scheme P there is a renormalized trajectory (RT), which 
starts at the ultra-violet flxed point. More than ten years ago, Iwasaki estimated a RT by 
matching Wilson loops (obtained in a certain approximation), and proposed to use it as 
an improved action Recently Hasenfratz and Niedermayer have reminded us that any 
point on the RT is a "perfect action" [^, since there the long distance behavior must be 
same as at the fixed point. Therefore if we find a RT, we have a perfect action. Even if it 
is an approximate one, it may serve as an improved action. 

In this paper we will determine the RG flow in the following approximation; the coupling 
space is restricted to two parameters and the action is assumed to have the following form, 

S = /?{cixi5](l-iReTrC/p,a9) 

plaq 

+ cix2^(l-^ReTrC/,ect)} 

red 

= -ia*(cixi + 8cix2)/3 5lTrF^-(^)' + 0(a')- (2) 

Here a is the lattice spacing and Upiaq and Urect correspond to 1 x 1 and 1x2 loops 
respectively. We impose the condition, cixi -I- 8cix2 = 1- 

2 Technique to determine the RG flow 

We adopt Swendsen's factor-two blocking scheme The blocked link variable is con- 

structed as: 

Qf,{x) = U^{x)U^{x + /,) + 1 ^ U,{x)U^{x -f v)U^{x + v + ^i)Ul{x + 2^l). (3) 

Projecting Q^{x) onto SU(3), we get a new blocked variable U' as max{i?e Tr Q^(x)?7^(x)}. 

From the original configurations, {?/}, generated by the action S with parameters 
(cixi,cix2) in Eq.(||), we obtain the blocked ones, {U'}, after blocking. If the blocked 
configurations, {C/'}, can be considered as generated by an action S' with parameters 
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(cixi',cix2') in Eq.(|), then (cixi,cix2) ^ (cixi', Cix2') can be regarded as the coupling 
flow associated with this blocking. 

For the determination of the effective action S' , we use a Schwinger-Dyson method 
This method is based on the following identity. For a link Ui, consider the quantity; 

< Im Tr(A"C/;Gn >= J DU Im Tii\''UiGf)e-^ /Z, (4) 

where A'^ stands for Gell-Mann matrices. G/ is a sum of "staples" Gj*^ for the link I, 
Gi = {Pj/6)G^ . The action S is assumed to have the form ^iRe TrlliGi. For the 
present analysis, 7 corresponds to a plaquette and a rectangle. Eq.(|4|) should be invariant 
under the change of variables [// ^ (1 + ieX'^)Ui. Setting the terms linear in e to zero, we 
get the identity, 

/ DU[ Re TriiX" fUiGf) + Im Tr(A"C/,GniniTr(A"C/,GO]e"'^ = 0. (5) 
Summing over c in the formula (^, we obtain the Schwinger-Dyson equation, 

^Rc<Tr(C/,Gf)>=^^ { -Rc < Tr(C/,Gf C/,G7) > (6) 
7 

+Re < Tr(Gr(G7)t) > - | < Im TT{UiGf)Im Tt{UiGI) >}. 

Here we used X^Li Tr(AM)Tr(A'=B) = 2TrAB ~ fTrATrB. We apply this equation to the 
blocked configurations, and calculate the expectation values < • • • > on both sides. Now 
Eq.(^ may be considered as a set of linear equations with /3^'s as unknowns. 

In general we may use other quantities instead of G", but the number of equations 
becomes equal to the number of unknowns if we take G" appearing in our action. In 
comparison with the demon method , larger loops such as Tr(G"(G7)^) are involved in 
the Schwinger-Dyson equation. 



3 Results of RG flow 

We have used lattices of size 8^ and 16^. About 2000 configurations separated by every 10 
sweeps are used to determine each parameter set (/3ixi) /3ix2) = /3(cixi, 01x2)- 

In Fig.|l|, the result of the coupling flow is shown. Arrows denote the measured coupling 
flow. We perform a Monte Carlo simulation with (cixi,cix2) (the starting point of an 
arrow), and then block the lattice. From the blocked links we estimate a blocked action 
(the end point) through the formula (^ above, by assuming that the blocked action has 
only two terms, plaquette (cixi') and rectangle (01x2')- It is remarkable that the arrows 
show a flow of simple structure. They flow into a trajectory. This is a strong evidence that 
the RT corresponding to this blocking scheme can be seen even in the two-coupling space. 

At the strong coupling region, the RT does not look like a linear curve but like a 
parabolic curve. This is completely consistent with the strong coupling expansion. In the 
strong coupling regime, (3 is small, and the Boltzmann weight e~^^ can be expanded as a 
polynomial in f3. If we include the 1x2 rectangle in the action in addition to the 1x1 
plaquette, then the first few loops can easily be evaluated in the regime where /3i x 1 , /9i x 2 <C 
1. In particular, the plaquette and the rectangle give: {Wn) « , and {W12) ~ 

(nr")^ + ^%t^- The string tension can then be evaluated at strong coupling as: 

2 , {W12) , /Axl /?1X2\ 
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Figure 1: Renormalization group flow for QCD in two coupling space 



Consider now the strong coupling limit /3ixi,/3ix2 ^ 0. The left-hand side should go to 
+00, indicating that the lattice spacing a diverges in the strong coupling limit |^3|. 

What is remarkable in the expression for the string tension is that it does not diverge 
unless + — > 0. In other words, as the renormalized trajectory approaches the 
high-temperature fixed point (HTFP) (/5ixi = /3ix2 = 0), it must approach the parabola 

/3ix2 = -^/3?xi- (8) 

As one goes away from the HTFP along the RT, the lattice spacing shrinks and larger 
loops must be considered to evaluate the string tension. A straightforward evaluation of the 
Creutz ratio X22 already shows that the parabola gets distorted "upwards" as one considers 
points on the RT further and furtheraway from the HTFP. 

On the other hand, in the vicinity of the Ultra- Violet Fixed Point (UVFP), where the 
lattice spacing is very small, the usual Symanzik improvement (tree-level, tadpole, one-loop 
or non-perturbative) should be effective at removing C(a^) cutoff effects. So one expects, 
as one approaches the UVFP, that the ratio — /3ix2//3ixi will approach the Symanzik value 
0.05. Our MCRG study shows that this ratio is approached from above, as seen also 
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when implementing the Symanzik scheme beyond tree-level. Here we exhibit the smooth 
connection from the UVFP to the HTFP in the two-coupling plane. 



4 Actions in two-coupling space 

The analysis in the previous section indicates that along the RT the ratio — Cix2/cixi 
should increase with the lattice spacing, at least until the lattice becomes very coarse. Here 
we study the rotational symmetry of the heavy quark potential of several approximations 
to the RT, characterized by a fixed value of cix2/cixi- 

We parameterize the two-dimensional plane by lines, cixi + 8cix2 = 1, with Cix2 
as a parameter; Cix2 — corresponds to the standard Wilson action, Cix2 = —1/12 to 
Symanzik's (t ree level with no tadpole improvement), Cix2 = — 0.331 to Iwasaki's H. 
We include one more steeper line with cix2 = —1.4088 which we call DBW2. This line 
goes through the point obtained after double blocking from a lattice with Wilson action 
iP = 6.3, 32^ X 64 ^ /3ixi ~ 7.986, /3ix2 ~ -0.9169) in 2-dimensional coupling space by 
the canonical demon method |l^ . 

First let us determine the lattice spacing. We measure the string tension, cTphysO,'^ 
from the heavy quark potential. The heavy quark potential is obtained from Wilson loops 
W{R,T) calculated from smeared links. The smearing technique proposed in ||l^ reduces 
contaminations from excited states to a very large extent. 

To determine the heavy quark potential V{R), links are first smeared along the spatial 
directions; W{R, T) is then computed in terms of the smeared links. The Wilson loops are 
expected to behave as 

W{R, T) = C{R) exp [-V{R)T] + ^ C,(i?) exp [-Vi{R)T] (9) 

where C{R) is the overlap function and V{R) is the heavy quark potential, the second term 
is the contribution from excited states. We obtain V{R) by applying least square fitting to 
the first term in Eq.(H) as a function of T for each R. The fitting range of T is chosen for 
each R such that the effective mass plots in T become approximately constant neglecting 
excited states in the fitting. The number of smearing steps is determined so that the overlap 
function C{R) is closest to one for each R. 

Next we extract the string tension from V{R). We take the following Ansatz, 

V{R) = A+^ + aR (10) 
R 



The string tension, a ~ cfphysci^ , is the coefficient of the linear term in Eq.(lO). We fit 
the heavy quark potential to Eq. ( p^ in the range where Ra corresponds approximately to 
physical distances about OAfm ^ 1.5/m. Statistical errors are estimated by the jackknife 
method with bin size one. We set -s/CTphiys ~ A2QMeV and fix the scale a. Results are 
summarized in Table |l]. 

Simulations are performed on a lattice of size 12"^ x 24 . Thermalization is 5000 sweeps, 
while the interval between Wilson loop measurements is 500 sweeps. The number of con- 
figurations used is 100. 

If the action lies on the RT, long range quantities behave in the same way as in the 
continuum. Therefore, even if the lattice is coarse, the rotational symmetry is expected to 
be recovered. Now we investigate the symmetry breaking in the (cixi,cix2) plane. 

We construct a quantity which measures the amount of rotational symmetry breaking; 



E 



[V{R) - Vor,{R)? / V- 1 



^ . V{R)HV{Rf \ 6V{RY 
on-axis Vott-axis 



(11) 
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action 


f3 cr a[fm] A 


Wilson action 
Clxl = 1 
Clx2 = 


5.45 0.463(45) 0.320(16) 0.082(16) 
5.55 0.2840(46) 0.2504(20) 0.049(10) 
5.65 0.1936(72) 0.2067(38) 0.0254(20) 


Symanzik action (tree) 

Clxl = 1 - 8cix2 

C1X2 - -1/12 


3.70 0.4666(89) 0.3209(31) 0.056(13) 
3.90 0.2634(41) 0.2411(19) 0.0317(53) 
4.10 0.1377(48) 0.1743(30) 0.0123(12) 


Iwasaki action 
cix2 = -0.331 


1 on c\ fi'ifi'i{70\ f^'i7AQ.^o^\ n c^'ifi^ f fi^\ 
l.yU U.0o0o(^ ( Z) U.oi 4o(^Zl ) U.UoD 1(^00 ) 

2.00 0.465(13) 0.3204(45) 0.0070(45) 
2.10 0.3496(79) 0.2778(31) 0.0114(27) 
2.20 0.2031(75) 0.2117(39) 0.0057(28) 
2.30 0.1333(63) 0.1715(41) 0.0071(28) 


DBW2 action 
C1X2 = -1.4088 


0.62 0.687(34) 0.3894(96) 0.0147(17) 
0.65 0.512(18) 0.3362(59) 0.0224(50) 
0.70 0.362(13) 0.2827(51) 0.0062(24) 
0.75 0.2489(71) 0.2344(33) 0.0054(16) 
0.80 0.1631(55) 0.1897(12) 0.0078(30) 



Table 1: String tensions and the lattice spacing determined from them 



where Von{R) is a fit to Eq.(pl]|) determined from on- axis data only. Here "on- axis data" are 
those on points along a lattice axis. dV{R) is the statistical error of V{R). We employ A as a 
measure of the difference between on-axis and off-axis potential, and therefore a measure of 
the rotational symmetry breaking. We calculate this quantity with three improved actions 
at various lattice spacings, i.e., various /3's. 

Results are shown in Table |l| and Fig.^. For small Cix2, like Wilson or tree-level 
Symanzik actions with no tadpole improvement, A is proportional to a^. As |cix2| in- 
creases, the rotational symmetry is better recovered. To the present accuracy and in the 
range of couplings studied, the breaking of the rotational symmetry can be already ignored 
at the level of Iwasaki's action. 

5 Concluding remarks 

In this paper, we calculate numerically the RG flow in two-couphng space (cixi,cix2), as 
the first step to study the RG flow. Although it is very remarkable that the flow structure 
can be already observed, the results may have "truncation" errors. After a blocking trans- 
formation, the action has in general infinitely many couplings. Therefore end points of the 
arrows in Fig.|^ generally deviate from the two-dimensional plane. To estimate this effect, 
we calculate the behavior of arrows, whose starting points are near the RT on the two- 
dimensional plane, in three-dimensional space: {(3ixi, f3ix2, (3twist), see FigJ|. Deviations 
from the (/3ixi, /?ix2) plane are small especially near the RT. This is expected also from 
the strong coupling expansion: because it takes 3 plaquettes to cover the "twist" loop, its 
effective coupling will behave like Pfy^i near the HTFP. 

It is very encouraging that here we have seen a clear flow behavior which strongly 
indicates the existence of a renormalized trajectory in the vicinity of the 2-coupling plane. 
We now plan to continue the Monte Carlo RG calculation (i) at higher /3's, in order to obtain 
the non-perturbative QCD P function, and (ii) in three- or four-coupling space including 
the chair and twist operators, to estimate truncation effects and evaluate the reliability of 
the results obtained here. 
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Figure 2: Rotational symmetry breaking 



All simulations have been done on CRAY J90 at Information Processing Center, Hi- 
roshima University, SX-4 at RCNP, Osaka university and on VPP500 at KEK (National 
Laboratory for High Energy Physics). H.M. would like to thank the Japan Society for the 
Promotion of Science for Young Scientists for financial support. 
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